Evidence for non-conventional pairing in Na^CoCVyH^O 
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We report on first investigation of the lattice dynamics in the novel superconducting material 
Nao.35CoO2 l.3H2O and the non-hydrated parent compound Nao.7Co02 by inelastic x-ray scatter- 
ing. The measured phonon dispersion along the F — M direction show a marked softening with 
hole doping of two optical phonon branches close to the Brillouin zone boundary. The phonon spec- 
tra, dispersion, and softening are well reproduced by first-principle calculations. The calculations 
indicates that the soft branches are mainly composed of Co-vibration modes. The estimation of 
the critical temperature based on electron-phonon coupling mechanism undisputedly points to a 
non-conventional superconducting state in this material. 
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As new superconductors are discovered, the identifica- 
tion of underlying pairing mechanisms is a key step in 
the understanding of the physics of these materials 
In the hydrated bilayer cobaltate Na^CoC^-yHbO su- 
perconductivity is observed in a narrow range of Na con- 
centration with an optimal T c of ^5 K The novel 
material is prepared from the hydration of the parent 
compound Na . 7 CoO 2 (NCO). The structure [HEll of 
Nao.7Co02 (space group P63/mmc) consists of alternat- 
ing layers of Na and C0O2 planes, stacked along the c 
axis in the hexagonal structure, showing some similarities 
with high-T c cuprates. The Na atoms can occupy two 
crystallographic sites (2b or 2d) , not necessarily in an or- 
dered pattern. In contrast to NCO, Nao.35CoO2T.3H2O 
(NCOH) has a considerably expanded c axis due to the 
intercalation of ice-like sheets in the structure. 

The parent unhydrated compounds Na x Co02 have a 
rich phase diagram, ranging from a Curie- Weiss metal 
(0.5 < x < 1), a charge transfer insulator (a; = 0.5) to 
a normal paramagnetic metal (x < 0.5) As it turns 
out, the Co electronic properties are not fully character- 
ized. In NCO and NCOH, Co supposedly coexists in two 
formal valent states, Co 3+ (5=0) and Co+ 4 (5=1/2), 
both in low-spin configuration. Thus, the Co ions form 
a frustrated two-dimensional triangular lattice as oppo- 
site to the Cu square lattice in cuprates. The overall 
picture is further complicated by the highly controver- 
sial electronic structure. A crucial issue is the existence 
of "hole pockets" on the Fermi surface : Hole doping 
results in the depletion of the Co-^2g conduction band. 
Under trigonal distortion, the t2 g triplet splits into two 
e' g bands and one a g band, which cross at the Fermi en- 
ergy. The latter forms a large hole-type Fermi surface 
around T, but the main contribution to the hole den- 



sity of states at the Fermi energy is due to the Co-e^ 
holes, supposedly distributed in 6 small pockets. These 
were predicted by calculations lalSlllfJ and, in the case 
of non-conventional pairing, their presence could help 
to distinguish between different pairing symmetries [ll|. 
However they have not be detected in photoemission ex- 
periments \VA . More recently, the formation of hole 
pockets are found robust against dynamical Coulomb cor- 
relations while no such a Fermi-surface topology is 
reported in the strong coupling limit 0] . The proximity 
of charge, magnetic [lij and structural 01 instabilities 
and the frustrated Co-0 triangular lattice suggests that 
the novel cobaltate could have exotic superconducting 
states. No information (experimental or theoretical) ex- 
ists on phonon dispersion and electron-phonon coupling 
in this material and the question remains of whether a 
conventional electron-phonon mechanism would yield the 
measured T c . 

In this work, we provide concrete evidence for the ex- 
otic character of the superconducting state. We measure 
the phonon dispersion in NCO and NCOH single crys- 
tals by inelastic x-ray scattering (IXS). The phonon dis- 
persion and the electron-phonon coupling are calculated 
from density functional theory. The results reveal pro- 
nounced softening of two optical phonon modes mainly 
formed by Co vibrations, while the calculated electron- 
phonon coupling strength definitively excludes conven- 
tional pairing mechanism. Further, the softening corre- 
sponds closely to what one would expect from a nesting 
mechanism arising from the e' g hole-pocket Fermi surface. 
These findings go a long way in identifying the character 
of the superconducting state in NCOH 

The small sample size (400 x 200 x 50 ^tm 3 ) excludes 
neutron scattering while the focused x-ray beam makes 
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IXS the privileged technique for measuring phonons . 
The experiment was carried on the ID-28 beamline at 
the European Synchrotron Radiation Facility (ESRF). 
Details about the experimental setup has been described 
elsewhere 19] . The resolution of 3 meV was achieved 
by selecting the Si(9,9,9) order of reflection for both 
the monochromator and the 5 spherical analyzers, and 
owing to the extreme backscattering geometry of the 
IXS spectrometer. The monochromatic beam is fo- 
cused by a toroidal mirror into a spot of 250 x 90 /im 2 
(horizontal x vertical) full width at half maximum at the 
sample position. NCO was mounted in a vacuum cham- 
ber at room temperature in order to reduced background 
scattering, and remove humidity from the air. To pre- 
serve hydration, NCOH was mounted in a closed-cycle 
cryostat and maintained at 100 K during the measure- 
ment. The lattice parameters for both samples were es- 
timated by monitoring the (1,0,0) and (1,0,1) Bragg re- 
flections on the detector. We obtained a = 2.82 A and 
c = 11.63 A in NCO and a = 2.83 A and c = 19.6 A 
in NCOH, in good agreement with the values reported 
in the literature for the two phases In the latter, we 
measured a critical temperature T c of 4.6 K by magnetic 
susceptibility. All the phonon spectra were measured in 
the second Brillouin zone around (1,0,0). To extract the 
phonon frequencies, the experimental spectra were least- 
square fitted to a sum of Lorenzian pairs, weighted by 
the Bose factor. 

Figure 1 shows the IXS phonon spectra in NCOH and 
NCO, along the T-M direction. Well defined phonon- 
like excitations are observed up to 80 meV, while non- 
dispersive features show up in most of the spectra up to 
the 10 meV region. As discussed hereafter, these are not 
correlated to any calculated phonon modes of NCOH and 
NCO, and may be reasonably attributed to the presence 
of impurities or to Na disorder in the sample. The spectra 
are compared to first- prin ciple calculations. Electronic 
structure calculations |2(| are performed using Density 
Functional Theory (DFT) in the local density approxi- 
mation (LDA) and ultrasoft pseudopotentials 21]. The 
electronic wavefunctions and the charge density are ex- 
panded using a 35 and a 350 Ryd cutoff respectively. The 
dynamical matrices and the electron-phonon coupling are 
calculated using Density Functional Perturbation Theory 
in the linear response |20j . For the electronic integration 
in the phonon calculation we use a N q = 6x6x2 uni- 
form k-point mesh and and Hermite-Gaussian smearing 
of 0.05 Ryd. For the calculation of the electron-phonon 
coupling we use a finer Nk — 25 x 25 x 25 Monkhorst 
Pack mesh. The electron-phonon coupling for a phonon 
mode v with momentum q is derived from: 
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FIG. 1: IXS spectra in non-hydrated and hydrated com- 
pounds measured along the T-M direction (open circles). The 
measurements are compared to ab-initio calculations of the 
dynamical structure factor (solid lines). The calculated spec- 
tra are shown for Na occupying sites 2b and 2d. The main 
dispersive feature is pointed by arrows. 



the electron density of states at the Fermi level and ekn 
are the DFT energy bands. The electron-phonon matrix 
element is fl£ njk+qm = (kn | SV/ fai q „ |k + qm) / v /2cj q „ , 
where u q „ and u)(\v are the amplitude of the displace- 
ment and the frequency of the phonon v of wavevec- 
tor q respectively while V is the Kohn-Sham poten- 
tial. The average electron-phonon coupling then reads 
as A = jj- v ^qv • The doping of NaCo02 is modeled 
simulating a charged cell with overall number of elec- 
trons corresponding to Naa;Co02. Since the positions of 
Na atoms are not determined we perform different simu- 
lations for Na occupying the 2b or the 2d positions. Cal- 
culating the phonon dispersion of Nao.3sCo02-yH2 is 
intractable due to the large number of atoms in the cell. 
However it has been shown [2^ that the main effect of 
the intercalation of water in Naa;Co02 is the c— axis ex- 
pansion. Unfortunately calculations of Na.,;Co02 with 
an expanded lattice parameters produces very flat bands 
with two very steep unphysical Na bands crossing the 
Fermi level |22| . As a consequence, the phonons were 
calcualted in Nao.3sCo02 and are here compared with 
the experiment in Nao.35Co02-yH20. 

Despite the approximations involved in the calcula- 
tions, there is broad agreement with experiment as ob- 
served in Fig. 1. The measured phonon linewidths are 
systematically found larger (5-6 meV) than the experi- 
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mental resolution (3 meV) . The broadening effect applies 
to all modes, independent of the momentum transfer, 
while calculated phonon linewidths are resolution lim- 
ited. This points to an extrinsic origin presumably re- 
lated to the crystal mosaicity (~ 1°), besides other cu- 
mulative effects. For comparison the calculated spectra 
were convoluted with a 6 meV broadening. 

We focus more particularly on the optical phonon 
mode (indicated by arrows in Fig. 1) which is dispersive 
along the T-M direction. The theoretical calculations 
correctly reproduce the frequency and spectral weight of 
this mode, although the calculated frequencies in NCO 
are slightly lower than measured ones close to M . This 
discrepancy can be easily understood when taking into 
account uncertainty in the determination of the Na con- 
centration (x — 0.7±0.1) and the strong sensitivity of the 
calculated phonon frequency of this mode on doping, in 
the vicinity of the M point. Using x = 0.8 for instance, 
we obtain u> = 47.6 meV at M, within 5% of the experi- 
ment. The same mode is seen in the hydrated sample (see 
arrows), but now significantly softened. At the M point, 
the experimental softening amounts to Aw ~ 7.6±1 meV. 

The experimental phonon frequencies are reported in 
Fig. 2 for NCO and NCOH along with the full calculated 
dispersion in model Nao.7Co02 and Nao.3sCo02 com- 
pounds. The phonon dispersion is obtained by Fourier 
interpolation of the dynamical matrices computed on the 
N q points mesh. The intensity of the dynamical struc- 
ture in the measured Brillouin zone is gray scale-coded. 
Raman and infrared data are also indicated as refer- 
ence. The theoretical phonon dispersions resemble each 
other closely except for the optical branches in the 40-45 
meV energy region which are sensibly softened on reduc- 
ing x from 0.7 to 0.35. The other branches are barely 
affected. The softening is consistently reproduced by the 
experimental results in the hydrated sample despite the 
mediocre quality of the spectra. Our calculation indi- 
cates that the softening is due principally to the reduc- 
tion of Na content and not to hydration. In addition 
to the calculations in Nao.3sCo02 with a reduced c pa- 
rameter (shown in Fig. 2a), the effect of hydration was 
modeled by simulating two layers of C0O2 separated by 
an expanded c parameter and having a formal charge 
corresponding to x = 0.35. The calculated phonon dis- 
persion (not shown here) exhibit a comparable softening 
of the optical branches, which signifies that softening is 
solely related to the doping in the C0O2 layers. As in- 
ferred from the eigenvector analysis of the computed dy- 
namical matrix, in the two compounds the high energy 
region of the phonon dispersion is mainly composed of 
oxygen vibrations while the intermediate energy region 
which shows the main softening effect is mainly due to 
cobalt vibrations. 

We can now estimate the critical temperature for 
Nao.7Co02 and for Nao.3sCo02 from the average 
electron-phonon coupling A and phonon frequency log- 




FIG. 2: Measured (open symbols) and calculated (solid cir- 
cles and dotted lines) phonon dispersion along high-symmetry 
directions in NCOH and NCO. For clarity purpose, only cal- 
culations with Na-26 sites occupied are shown. The structure 
factor intensity is indicated in normalized unit by the gray 
scale form black (1) to light gray (10 -6 ), and below 10~ 6 as 
dotted lines. Raman and infrared data are borrowed from 
Ref. I2H Left panels represent the projected phonon density 
of states. 



arithmic average (uj). The calculated values are A = 0.21 
and (ui) = 42.5 meV for x = 0.7 and A = 0.29 and 
(lu) = 40.3 meV for x = 0.35. Using McMillan's formula 
we obtain T c » 0.001 K for x = 0.7 and T c = 0.15 K 
for x — 0.35 with the Coulomb pseudo-potential con- 
stant set at fi* = 0.1. These are upper estimates of T c 
since in both systems /j,* is probably larger. Calcula- 
tions for the expanded C0O2 bilayer lead to A = 0.28, 
(u>) = 38.6 meV and T c = 0.12 K, showing that the ex- 
pansion of the lattice spacing weakly affects the electron- 
phonon coupling. Some earlier works indicate that water 
intercalation affects the valence of Co @, 0] ■ The agree- 
ment between the phonon spectra and our calculation 
indicates that this effect if any is negligible. Since the 
calculated T c , based on a conventional electron-phonon 
mechanism, gives results which are more than an order 
of magnitude smaller than the measured T c , we have a 
strong argument for invoking a non-conventional pairing 
mechanism. As noted in Ref. 111! the determination of 



4 



the pairing symmetry also requires a knowledge of the 
electronic structure. In particular a crucial aspect is the 
presence of hole pockets on the Fermi surface. As it turns 
out, our phonon spectra reveal information about these 
hole pockets through their influence on phonon frequen- 
cies. Indeed the softening of the optical phonon-mode is 
connected to the real part of the phonon self-energy due 
to the electron-phonon interaction, namely: 
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where /km stands for the Fermi function. If the 
matrix elements gfc n k+qm are assumed to be constant, 
IT„(q, ujqv) is related the one-electron susceptibility as 
calculated in Ref. 0- For x — 0.35 the real part of 
n„(q, ujqu) displays pronounced features at the M point 
which is close to the hole pocket nesting vectors. The 
net effect on the phonon spectrum is the softening of the 
phonon branches mentioned earlier. 

In conclusion, we present first measurements and cal- 
culations of the phonon dispersion of Nao.7Co02 and 
Nao.asCoCVyEbO using inelastic scattering of hard x- 
rays as a probe of bulk properties. A good agreement 
is found between theoretical and experimental phonon 
dispersion. The estimation of the superconducting criti- 
cal temperature assuming a phonon-mediated mechanism 
leads to values which are more than an order of mag- 
nitude smaller than the measured critical temperature. 
Moreover the measured and calculated softening of opti- 
cal phonon modes can be traced back to the existence of 
the hole pocket Fermi surface. These results point in the 
direction of a non-conventional superconducting state in 
Na .35CoO 2 -yH 2 O. 
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